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RESUMEN 
ABSTRACT 

We used galaxy evolutionary models in a hierarchical inside-out disk forma- 
tion scenario to study the origin of the main local and global properties of disk 
galaxies as well as some of the correlations between their global properties. We 
found that most of disk galaxy properties and their correlations are the result of 
three (cosmological) initial factors of the system and their dispersions: the virial 
mass, the halo mass aggregation history (MAH), and the angular momentum given 
through the spin parameter A. The MAH determines mainly the dark halo structure 
and the integral color indexes while A determines mainly the disk surface bright- 
ness and the bulge-to-disk ratio. We calculated star formation using a gravitational 
instability criterion and a self-regulation mechanism driven by the energy balance 
in the disk turbulent interstellar medium. The efficiency of star formation in this 
model is almost independent from the mass of the system. Therefore, the galaxy 
integral color index B — V oi the stellar mass-to-luminosity ratio, M/Lb, do not 
depend upon mass. We show that the luminosity-dependent dust absorption em- 
pirically determined by Wang & Heckman explains the observed color-magnitude 
and color TuUy-Fisher (TF) relations without the necessity of introducing a mass- 
dependent star formation efficiency. The star formation history in our models is 
driven by both the MAH and the surface density (A). The disks in centrifugal 
equilibrium form within growing cold dark matter halos with a gas accretion rate 
proportional to the halo MAH. The disks present exponential surface density and 
brightness profiles, negative radial color index gradients, and nearly flat rotation 
curves. We also calculated the secular formation of a bulge due to gravitational 
instabilities in the stellar disk. The intensive properties of our models agree very 
well with the observational data and they follow the same trends galaxies follow 
across the Hubble sequence. The infrared TF relations and the luminosity-radius 
relation predicted for the cosmological model used here also agree with the observa- 
tional data. The main shortcomings of the inside-out hierarchical formation models 
presented in this paper are the excessive radial color gradients and the dark halo 
dominion in the rotation curve decompositions. 
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1. INTRODUCTION 

The study of galaxy formation and evolution may 
be approached in a deductive or in an inductive way 
(e.g., Renzini 1994). In the former case, starting 
from some ab-initio conditions given by a theory of 
cosmic structure formation, one tries to model the 
galaxy evolution in order to predict the observable 
properties of the galaxies. In the latter case, start- 
ing from the present-day properties of galaxies, and 
through galaxy evolutionary models, one tries to re- 
construct the initial conditions of galaxy formation. 
In this paper we introduce models of disk galaxy for- 
mation and evolution based on both approaches. 

Deductive approach. Most of current mod- 
els about cosmic structure formation are based on 
the gravitational paradigm and the inflationary cold 
dark matter (CDM) cosmology. In these models, cos- 
mic structures build up hierarchically by a continu- 
ous mass aggregation process. Prom the point of 
view of galaxy cosmogony, a key question is whether 
this aggregation occurs through violent mergers of 
collapsed substructures or through a gentle process 
of mass accretion. According to this, two could be 
the most general scenarios of large-galaxy forma- 
tion, (i) In one case, the main properties of luminous 
galaxies, including those which define their morpho- 
logical types, are supposed to be basically result of 
a given sequence of mergers. This picture — that we 
shall call the rn.erger scenario — has been widely ap- 
plied in semianalytical models of galaxy formation, 
where the global properties of galaxies are calculated 
from the cosmological dark matter (DM) halo merg- 
ing histories using preconceived recipes for the gas 
cooling, disk mergers, star formation (SF), super- 
nova (SN) feedback, etc. (e.g., Lacey et al. 1993; 
Kauffmann, White, & Guiderdoni 1993; Cole et al. 
1994; Heyl et al. 1995; Kauffmann 1995; 1996; 
Baugh, Cole, & Frenk 1996; Baugh et al. 1997; 
Somerville & Primack 1999). (ii) In the other case, 
the formation and evolution of galaxies is related 
to a gentle and coherent process of mass aggrega- 
tion dictated by the shapes of the density profiles 
of the primordial fluctuations: galaxies coherently 
grow insidc-out. Wc shall call this picture — flrstly 
developed by Gmm (1981, 1987) and by Ryden & 
Gunn (1987) the extended collapse scenario. Since 
disk galaxies (^ 80% of present-day normal galaxies) 
could not have suffered major mergers due to the dy- 
namical fragility of the disks (Toth & Ostriker 1992; 
Weinberg 1998), the extended collapse scenario ap- 
pears to be more appropriate to study their evolu- 
tion. 

The question of whether the major mergers dom- 



inate or not in the assembly of protogalaxies depends 
on the statistical properties of the density fluctuation 
fleld. However, even if the DM structures assem- 
ble through chaotic and violent mergers of substruc- 
tures, the baryon gas tends to aggregate and collapse 
in a more (spatially) uniform fashion than DM do it 
(e.g., Kepncr 1997). On the other hand, the rcioniza- 
tion of the universe probably heated the gas at high 
redshifts to temperatures enough to prevent further 
its collapse within the small DM substructures of the 
hierarchy (e.g., Blanchard, Valls-Gabaud, & Mamon 
1992; Prunet & Blanchard 1999). This way, large 
galaxies could have formed smoothly with the gas 
trapped into the potential wells of their surrounding 
DM halos, avoiding an assembly through the merg- 
ing of baryon substructures. Only in high density 
regions of the universe, as in clusters of galaxies, 
luminous galaxies could have efficiently interacted 
among each other, thickening and disrupting their 
disks and probably forming lenticular and elliptical 
systems. Disk galaxies indeed typically are isolated 
objects. On the other hand, the interpretation of 
high velocity clouds as build-blocks of the formation 
of galaxies in the Local Group (Blitz et al. 1999) is 
one possible observational evidence for the extended 
collapse scenario: these clouds might constitute large 
reservoirs of cold gas which is smoothly accreted onto 
the galaxies (Lopez-Corredoira, Beckman, & Casuso 
1999). 

Inductive approach. Galaxy evolutionary 
models have shown that due to the relatively rapid 

disk gas consumption in SF, closed models arc not 
able to explain several properties of disk galaxies, 
as well as the wide range of colors, gas fractions, 
etc. that galaxies present across the Hubble se- 
quence (e.g., Larson & Tinsley 1978; Tinsley 1980; 
Larson, Tinsley, & Caldwell 1980; Kennicutt 1983; 
Gallagher, Hunter, & Tutukov 1984; Firmani & Tu- 
tukov 1992, 1994). On the other hand, it was shown 
that the SF time scale in disk galaxies is not mainly 
controlled by the initial gas surface density (Ken- 
nicutt 1983; Firmani & Tutukov 1994; Kennicutt, 
Tamblyn, & Congdon 1994). Hence, models where 
accretion of fresh gas is introduced are more realis- 
tic. What is more, in order to reproduce the radial 
abundance gradient of the Galaxy, the gas accretion 
should be such as an inside-out disk formation sce- 
nario (e.g., Portinari & Chiosi 1999, and more ref- 
erences therein). Gas accretion might also be nec- 
essary to maintain spiral structures. In the case of 
open (infall) models, galaxy formation and evolution 
could be two related processes where the SF time 
scale is mainly driven by the gas accretion rate at 
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which the disk is being built up. 

Gas infall and SF. The gas infall rate onto 
the disks may bo controlled by the global process 
of galaxy formation (cosmological accretion) and/or 
by a self-regulated process of SF. This latter pro- 
cess, proposed in the seminal papers by White & 
Rees (1978) and White & Prenk (1991), is commonly 
applied in the merger scenario models. According 
to this mechanism, the gas accretion rate is driven 
by the cooling of the hot gas corona sustained by 
the SN-injected energy. However, the self-regulated 
disk-halo SF models suffer from some inconsisten- 
cies. As Nulsen & Fabian (1996) pointed out, SN 
feedback over large scales occurs on roughly the same 
time scales as the SF, not fast enough to tightly regu- 
late the SF rate. On the other hand, because the disk 
turbulent interstellar medium is a very dissipative 
system (Avila- Reese & Vazquez-Semadeni 2000), gas 
and energy outflows tend to be confined into a re- 
gion close to the disk. The thick gaseous disk and 
the global magnetic fields are efficient shields that 
prevent any feedback toward the halo on large scales 
(e.g., Mac Low, McCray, & Norman 1989; Slavin & 
Cox 1992; Franco et al. 1995). The lack of observa- 
tional evidence for significant amounts of hot gas in 
the halos of disk galaxies confirms this picture and 
suggests that the halo-disk connection is not enough 
to self-regulate and drive the disk SF at the halo 
scales. 

An unified scenario of galaxy formation. 
The galaxy infall models suggested by the inductive 

approach can in a natural fashion be incorporated 
into the cosmological (deductive) extended collapse 
scenario. In Avila-Reese, Firmani, & Hernandez 
(1998, hereafter AFH98), starting from the cos- 
mological settings, the formation and evolution of 
galaxy DM halos were calculated through a sim- 
plified approach whose results are in good agree- 
ment with those of the cosmological N-body simula- 
tions (see also Avila-Reese et al. 1999). In Firmani 
& Avila-Reese (2000a, hereafter FA), the formation 
and evolution of baryon disks within the growing ha- 
los produced with this approach was calculated in 
order to study the structural and dynamical proper- 
ties of disk galaxies. Here, we apply a similar method 
including bulge formation, with the aim to explore 
whether the origin of the main properties of disk 
galaxies and their correlations (the Hubble sequence, 
for example) can be related or not to cosmological 
initial conditions and factors. In an acompanying 
paper (Firmani & Avila-Reese 2000b) results about 
the evolution of the model galaxies will be presented 
and compared with the observational data. 



In the extended collapse scenario we develop 
here, the relevant assumptions we introduce are: (1) 
spherical symmetry and adiabatic invariance in the 
gravitational collapse of DM halos, (2) spin param- 
eter A produced by cosmological torques and con- 
stant in time for a given halo, (3) gentle incorpora- 
tion of the gas (no mergers) into a disk in centrifugal 
equilibrium, (4) detailed angular momentum conser- 
vation and adiabatic invariance during the gas col- 
lapse, (5) stationary SF induced by global disk gravi- 
tational instabilities and regulated by an energy bal- 
ance within the disk interstellar medium, and (6) for- 
mation of a bulge due to gravitational instabilities in 
the stellar disk (secular formation scenario). 

A consequence of item (5) is that in our mod- 
els the SF efficiency is almost independent from the 
galaxy mass or circular velocity. An important ques- 
tion is whether the observed "color" Tully-Fisher or 
color-magnitude relations can be explained without 
invoking SF efficiency and/or disk-to-total mass ra- 
tios depending upon mass (luminosity) . We propose 
that the observed luminosity-dependent dust absorp- 
tion (e.g., Wang & Heckman 1996) explain these em- 
pirical relations. 

In the last years, several methods to model galaxy 
formation and evolution were proposed: the numeri- 
cal simulations (N-body-|- hydrodynamics), the semi- 
analytical method, and the analytical method. We 
have called our method seminumerical, essentially 
with the aim to differentiate it from the other ones. 
In fact, the virialization of the halos is calculated 
with an itterativc procedure and the disk formation 
and evolution is followed by solving the hydrody- 
namical equations (in the radial direction) with a full 
numerical method. As in the numerical simulations, 
with the seminumerical method, we try to model the 
evolution of the internal structural, dynamical, kine- 
matical and SF processes of an individual galaxy 
(dark-|- luminous). On the other hand, we are also 
interested in calculating statistical properties of a 
whole population of galaxies. The techniques used in 
semianalytical and analytical methods have proved 
to be especially powerful in this undertaking. In our 
method we use some related techniques. It is im- 
portant to mention that, from the physical point of 
view, our seminumerical method and the analytical 
method are similar: both are based on the extended 
collapse scenario. The results obtained with the ana- 
lytical method in recent papers (Dalcanton, Spergel, 
& Summers 1997; Mo, Mao, & White 1998; van den 
Bosch 1998, 1999) encourage us to study in more 
detail this scenario, following the overall process of 
formation and evolution of individual disk galaxies 
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(dark+himinous) . 

The plan of the paper is as follows: in § 2 we 
describe our method and explain the strategy in or- 
der to calculate the galaxy models; in §§ 3 and 4 we 
present the local and global properties of the model 
galaxies at 2; = 0, respectively, and compare them 
with the observational data; in § 5 wc stiidy some 
important correlations among the global properties 
of the models, particularly those which go across the 
Hublc sequence, and wc explore whether the "color" 
TuUy-Fisher and color-magnitude relations can be 
explained by the luminosity-dependent dust absort- 
pion; the fundamental factors of the galaxy models 
and the way they determine the galaxy properties 
are presented in § 6; we discuss the main difficulties 
found in the extended collapse scenario in § 7, and 
finally, in § 8, we summarize the paper and present 
the conclusions. 

2. THE METHOD 

Following, we briefly describe the seminumerical 
method used to calculate the formation and evolu- 
tion of isolated disk galaxies in the extended col- 
lapse scenario based on the hierarchical formation 
picture. Most of the steps of this method were al- 
ready described in previous works (Firmani et al. 
1996; Avila-Reese 1998; AFH98; FA). 

Virialization of the DM halos. The first step 
in our model is to calculate the gravitational collapse 
and virialization of isolated DM halos starting from 
a primordial density fluctuation field. This was done 
in AFH98 and in Avila-Recsc 1998. In order to gen- 
erate the hierarchical MAHs of the DM halos we use 
the extended Press-Schechter approximation based 
on the conditional probabilities for a Gaussian ran- 
dom fleld (Bower 1991; Bond et al. 1991; Lacey & 
Cole 1993). For a given present-day mass Mnom, we 
generate a set of MAHs through Monte Carlo simula- 
tions. We follow the aggregation history of the main 
progenitor by identifying the most massive subunit 
of the distribution at each time. The gravitational 
collapse and virialization of the DM halos is calcu- 
lated with an iterative seminumerical method, as- 
suming spherical symmetry and adiabatic invariance 
during the collapse. This method is based on the 
secondary infall model (e.g., Gunn 1977; Zaroubi & 
Hoffman 1993), and it expands upon it by allowing 
non-radial motions and arbitrary initial conditions 
(MAHs in our case). 

Not all the mass Mnom virializes at z = 0. The 
mass shells that are already virialized roughly corre- 
spond to those within the virial radius at which 
the halo mean overdensity drops below the critical 



value Ac given by the spherical collapse model. For 
the cosmological model used here Ac = 179 (e.g., 
Padmanabhan 1993). Analysis of halos identified 
in numerical simulations confirm this result (Cole & 
Lacey 1996; Eke, Navarro, & Frenk 1998). The mass 
contained within r„ is the virial mass My which, de- 
pending upon the MAH, is equal to 0.7-0.9 times 
Mnom (see also Kufl 1999). 

The halos calculated with our method present a 
diversity of density profiles, the most typical being 
close to that suggested by Navarro, Frenk, &: White 
(1996, 1997) (AFH98). This diversity in the struc- 
ture of the halos is related to the dispersion in the 
MAHs. As it was tested for a flat CDM model with 
cosmological constant, the diverse density profiles 
obtained with our method agree rather well with the 
results of cosmological N-body simulations (Avila- 
Reese et al. 1999). 

Disk build-up. Baryon matter is able to cool 
and further collapse within the virialized DM halo. 
The collapse can be halted only by centrifugal forces 
or when gas is transformed into stars. While the sec- 
ond case could be related to the origin of large bulges 
and elliptical galaxies, the former is commonly in- 
voked to explain the formation of disk galaxies. The 
build-up of disks within the evolving DM halos is 
carried out as follows: 

(1) We consider that, at the beginning, baryon 
matter has mass and angular momentum distribu- 
tions similar to those of the DM; it is assumed that 
the spherical shells are in solid body rotation. 

(2) Once the current mass shell has attained its 
maximum expansion radius, a fraction fa of its mass 
is transferred in a virialization time into a disk in 
centrifugal equilibrium. This mass fraction is as- 
sumed to be only in form of gas. Since for galaxy 
halos the time scale of gas cooling is generally smaller 
than the dynamical time scale, the shell virialization 
time is a reasonable — probably minimum — time 
scale for gas incorporation into the disk, considering 
that the gas will be shock-heated during the virial- 
ization process. The radial mass distribution of the 
disk gas is calculated equating its specific angular 
momentum to the one of its final circular orbit (de- 
tailed angular momentum conservation is assumed). 
The specific angular momentum jsh acqiiired by each 
collapsing shell during the linear regime is estimated 
under the assumption of a constant spin parameter 
A (= {J\E\^^^)/{GM^/^)): 

dJjt) _ GM{tf/^ 1 d\E{t)W 
dM{t) \E{t)\^/^y'2M{t) 2dM{t)) 

(1) 
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where M , J and E arc the total mass, angular mo- 
mentum and energy of the object at time t, respec- 
tively. According to analytical studies and to the 
outcomes of cosmological N-body simulations, the 
DM halos have values of A given by a log normal 
distribution with an average of ~0.05 and a width in 
the logarithm less than one (e.g., Catelan & Theuns 
1996 and references therein). 

(3) The gravitational drag on the total system 
produced by the central disk is calculated at each 
time with the adiabatic invariant formalism (e.g., 
Flores et al. 1993). 

Disk star formation and hydrodynamics. 
In our model, the SF and the internal hydrodynamics 
of the disk are regulated by an energy balance be- 
tween the SN and gas accretion kinetic energy input, 
and the turbulent energy gas dissipation (Firmani et 
al. 1996). Star formation is turned on at radius r 
when the Toomre gravitational instability parameter 
for the gas disk. 



falls below a given threshold (Toomre 1964); k is the 
epicyclic frequency, Vg the gas rms turbulent veloc- 
ity, and Tig the gas surface density. Thus, the SF is 
controlled by a feedback mechanism such that, when 
a gas disk column is overheated by the SF activity, 
SF is inhibited and the disk column dissipates the ex- 
cess energy to lower Vg back to the value determined 
from the Toomre criterion threshold. Owing to sev- 
eral theoretical and observational pieces of evidence 
(see references in FA) we fix Qg = 2 instead of the 
value of 1 given by Toomre (1964) for an infinitely 
thin disk. While the SF rate is rather insensitive 
to Qg, the disk thickness strongly depends upon it. 
When Qg = 2 is used, for a model of the Galaxy 
we obtain the observed gas and stellar thicknesses. 
The gas loss from stars is also included. The gravi- 
tational dynamics of the evolving star and gas disks, 
and the DM halo are treated in detail. 

The local disk galaxy models presented in Fir- 
mani & Tutukov (1994) included the integration of 
luminosities in B and V bands for all the stellar pop- 
ulation. A Salpeter initial mass function with a min- 
imal star mass of O.IMq, and solar metallicities were 
used. Here, the results from these models are used 
to calculate the surface B— and y— band brightness 
at every radius and throughout the evolution of the 
galaxy models (see Firmani et al. 1996). We find 
that the approximation to the population synthesis 
we use provides B — V colors that can be less by 
~ 0.1 mag in the red, and by ~ 0.05 mag in the 



blue than the respective values given by more so- 
phisticated models (e.g., Bruzual & Chariot 1993; 
Chariot, Worthey, & Bressan 1996). 

According to our scheme, time by time and at 
each radius, the growing disk is characterized by the 
infall rate of fresh gas by unit of area, Eg(r, t), the 
gas and stellar disk surface density profiles, Eg(r, t) 
and T,s{r,t), the total rotation curve (including the 
growing DM halo component), Vr{r,t), and the SF 
rate (r, t) determined by the energy balance in 
the vertical gaseous disk and by a Toomre criterion. 
It is important to note that in our self-regulating 
SF mechanism the feedback happens only within the 
disk and not at the level of the whole halo. This 
is justifiable, as was already explained in the Intro- 
duction, because the galaxy disks are dense and very 
dissipative systems that confine gas and energy out- 
flows. 

Bulge formation. Recent observational and 
theoretical studies are challenging common precon- 
ceptions about galactic bulges; for a review see 
Wyse, Gilmore, & Franx (1997). These studies tend 
to show that several bulge formation mechanisms 
could be working in galaxies. Bulges in galaxies with 
low and intermediate bulge-to-disk (b/d) ratios may 
be formed through secular dynamical evolution of 
the disks, whereas bulges in galaxies with large b/d 
ratios could have been formed separately from the 
disks, through an early dissipative collapse, and/or 
from mergers. 

In order to estimate the bulge mass we have in- 
troduced a recipe according to which the stars of the 
central "cold" disk region where the stellar Toomre 
instability parameter Qs(r) = ■/36Gs^(r) ^^^^ than 
1, are transferred to a spherical component in such 
a way that Qs remains equal to 1; Vg and Eg are the 
radial stellar velocity dispersion and stellar surface 
density. The physical sense of this recipe is in agree- 
ment with the secular scenario of bulge formation 
where gravitational instabilities in the stellar disks 
produce bars which dissolve forming a "hot" com- 
ponent, the bulge (see Norman, Sellwood, & Hasan 
1996, and references therein). Indeed, the similar- 
ity found in colors (Peletier & Balcells 1996), and 
scalelengths (de Jong 1996a; Courteau, de Jong, & 
Broeils 1997) between disk and bulge for late-type 
galaxies might mean that the disk and the bulge are 
closely related. It is worth emphasizing that our 
recipe of bulge formation is a very crude approxi- 
mation to a complex phenomenon about which, in 
fact, not much is known, both from the theoretical 
and observational points of view. 
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2.1. Selection of the models 

The seminumerical method described above al- 
lows us to calculate the properties of disk galaxies 
at any time, particularly at the present-day, start- 
ing from the cosmological initial conditions. Since 
our aim is to study general behaviors, here we shall 
use only one representative cosmology, the standard 
CDM model with the following parameters: matter 
density Qm = 1, baryon fraction density flf, — 0.05, 
Hubble parameter Hq = h x 100 kms~^Mpc^^ with 
h = 0.5, the amplitude of fluctuations on 8/i~^Mpc 
scale as = 0.6. At galaxy scales, the power spectrum 
of this model is similar to that of the popular flat 
CDM model with cosmological constant (57a ~ 0.7, 
h « 0.7) normalized to the COBE— satellite mea- 
surements. We calculate the power spectrum ac- 
cording to Sugiyama (1996). Once the cosmologi- 
cal model is fixed, the galaxy model to be studied 
should be defined. A given model is characterized 
by its (dark-|-baryon) virialized mass M^Q, its MAH, 
its initial spin parameter A, and by the mass fraction 
fd of My that is incorporated into the disk. In order 
to obtain a complete galaxy population, one should 
calculate models with masses, MAHs, and A's drawn 
from the corresponding statistical distributions given 
by the cosmological framework. A complete model- 
ing of the disk galaxy population, however, is beyond 
the scope of this paper. Regarding fd-, it is proba- 
bly the same for all the galaxies, although some as- 
trophysical process could change their initial values 
(e.g., Natarajan 1999). Here we shall assume fd con- 
stant and equal to 0.05 (see Mo et al. 1998 and FA 
for some justifications for taking this value). This 
value is within the range predicted by the theory of 
light elements production, ^bh^.^ ~ 0.025 — 0.053, 
as well as by the range given by the global budget of 
baryons, f^b/iQ.s « 0.014 - 0.080 (Fukugita, Hogan, 
& Peebles 1998). Note that, taking fd ^ 0.05 for the 
standard CDM model, we are assuming that in the 
galaxy halos almost all the baryon matter is incor- 
porated into the disk. 

The aim of this paper is to explore to what ex- 
tent the extended collapse scenario is able to pro- 
vide the correct initial conditions for galaxy evolu- 
tionary models. With this aim in mind, we shall 
study and compare with observations (i) the main 
structural and luminosity characteristics of galaxies 



^ As it was explained above, we give as initial parameter the 
nominal mass Mnom- Nevertheless, not all this mass is virial- 
ized; the most outer mass shells that encompase this mass are 
still in the process of virialization. We use the virial mass A/u 
as the current mass of the halo. This mass is approximately 
0.7 - 0.9A/nom, depending on the MAH. 



corresponding to a local disk galaxy population, and 
(ii) the main correlations among the global proper- 
ties, particularly those which go across the Hubble 
sequence. To achieve the former item, it is enough 
to calculate models for a significant range of masses, 
MAHs and A's, and test whether the obtained galaxy 
properties are realistic for these ranges. Regarding 
the latter item, as was mentioned above, we shall not 
model a complete galaxy population; only general 
trends in the correlations among the global galaxy 
properties will be obtained. As a matter of fact, 
the available observational galaxy samples are not 
complete enough so as to provide complete statisti- 
cal information on galaxy properties and their cor- 
relations. Only certain relations, as the TuUy-Fisher 
relation (TFR), were observationally studied with 
enough detail as to have a statistical description. 

We calculate models for three representative 
masses: M„om = SxlO^^Mg, 5xlO"M0, and 
5x10^^Mq. For each mass, three MAHs are se- 
lected from its distribution: the average, and two 
symmetrical deviations such that roughly 80% of all 
the MAHs are contained between them (see AFH98). 
For each mass and for each MAH, we calculate mod- 
els taking three values of A : 0.035, 0.05, and 0.1. 
Thus, we obtain a grid of 27 models. 

3. THE LOCAL PROPERTIES OF THE 
MODELS 

Local star formation rate. We obtain that 

SsFR(r) cc E^(r) (3) 

with n « 2, for most of the models and along a 
major portion of the disks, where Ssfr is the SF 
rate per unit of area (see also Firmani et al. 1996). 
This SF law is basically a consequence of the self- 
regulation feedback mechanism used in our models 
(Firmani & Avila-Reese, in preparation). Other self- 
regulation SF models also yield SF rates with (near 
to 2) power-law dependences on local Sg (see ref- 
erences in Struck & Smith 1999). Unfortunately, 
the observational data on the local SF rate of disks 
galaxies are still very poor. The data presented in 
Kennicutt (1989, 1998) show that the disks at inter- 
medium radii typically have a n « 2 SF law, while at 
inner radii the index is typically smaller than 2 and 
near to 1. Our SF prescription in inside-out forming 
disks is just a simple approximation to very complex 
phenomena; at this level, we think that our model 
predictions are succesful. 

Surface brightness profiles. The mechanism 
of disk formation used here naturally leads to an ex- 
ponential radial distribution of the disk surface den- 
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sity. Note that wc have not assumed initial angular 
momentuin distribution corresponding to an object 
in solid body rotation (e.g., Fall & Efstathiou 1980; 
Dalcanton et al. 1997). The crucial point is that we 
have assumed A constant in time (see eq. 1). Analy- 
sis of cosmological N-body simulations seem to con- 
firm this assumption (Gottlober, personal commimi- 
cation). In Figure la we show the surface brightness 
(SB) profiles in the B band for different models of 
A^nom =5xlO-^^M0. The SB profiles arc nearly ex- 
ponential and the models with large A form low SB 
galaxies (see also Dalcanton et al. 1997; Jimenez et 
al. 1998). For other masses the situation is the same, 
although with the mass the SB increases on average. 

Radial distribution of the B — V color in- 
dex. A natural consequence of the inside-out disk 
galaxy formation mechanism is the existence of neg- 
ative radial color index gradients along the disk. In 
Figure lb we plot the radial B — V color distribu- 
tions for the same models of Figure la. Multiwave- 
length surface photometry studies of galaxies (e.g., 
de Jong 1995,1996b) confirm that galaxy disks tend 
to be bluer at the periphery, although these gradi- 
ents are typically smaller than those predicted by 
our models. We shall return to this point in §6.2. 

Rotation curves. In Figure 2. the rotation 
curves for the same models of Figure 1 are depicted. 
The shape of the rotation curve depends mainly 
upon the spin parameter A and the MAH, and it 
correlates with the central SB. Observational studies 
seem to confirm the fact that galaxies with higher SB 
have more peaked rotation curves (e.g., Cascrtano & 
van Gorkom 1991; Tully & Verheijen 1997; Verhei- 
jen 1997). Models with the early active MAH have 
more concentrated rotation curves than models with 
an extended MAH. In Figure 2, the rotation curves 
of the Mnom = 5 X lOi°M0 and 5x10^^^^ j^od- 
els, corresponding to the average MAHs and A, are 
also depicted. At the Holmberg radius the less mas- 
sive models have nearly flat rotation curves, while 
the more massive galaxies present decreasing rota- 
tion curves. A similar trend was found from a large 
observational sample (Persic, Salucci, & Stel 1996). 
This could seem contradictory with the result that 
less massive DM halos are more concentrated than 
the more massive ones (Navarro et al. 1996,1997; 
AFH98; Avila-Reese et al. 1999). However, in our 
models, this is compensated by the fact that the disks 
in massive halos are more concentrated than in small 
halos (sec also Dalcanton et al. 1997), with the result 
that their contribution to the total rotation curve is 
more significant in the more massive halos than in 
the less massive ones. The disk mass fraction 



also influences on the rotation curve shape. If we fix 
fd < 0.05 then the rotation curves of all the model 
result less peaked than those presented in Figure 2 
(sec also Mo et al. 1998, FA). 

Regarding the decomposition of the rotation 
curves, as it was already shown in AFH98 and FA, 
we obtain that the halo component dominates almost 
everywhere for most of the models. This result is in 
potential disagreement with the rotation curve de- 
compositions inferred from observations (e.g., Carig- 
nan & Freeman 1985; Sancisi & van Albada 1987; 
Begeman 1987; Verheijen 1997; Corsini et al. 1998). 
In §6.3 we shall return to this question. 

4. THE GLOBAL PROPERTIES OF THE 

MODELS 

Among the global properties that our models pre- 
dict for a disk galaxy, we shall take into account for 
this study the integral B — V color index (calculated 
within a Holmberg radius), the i?-band luminosity 
Lb, the B-band exponential disk scale length hs, 
the B-band central SB ij,b„ or Eb„ (/Ub„ is given in 
magnitudes per arcsec^, and in Lbq per pc~^), 
the disk gas fraction fg (= ^ — ), the stellar 
bulge-to-disk ratio b/d, and the maximum rotation 
velocity Vm ■ The correlation matrix from a principal 
component analysis of these properties for the 27 
models calculated here is given in Table 1. The cor- 
relation coefficients give only a qualitative estimate 
of the correlations; note that the initial parameters 
of the models were not drawn randomly from the 
corresponding statistical distributions (§2.1). 

From Table 1 is seen that the observational 
trends across the Hubble sequence are reproduced 
for the intensive properties: the redder and more 
concentrated is the disk, the smaller is the gas frac- 
tion and the larger is the b/d ratio. Moreover, 
the models seem to populate a planar region in the 
HBo ~ {B — V) — fg or jiBo — {B — V) — b/d spaces, in 
agreement with observations (McGaugh & de Blok 
1997). B — V and fiBa ^^re almost independent one 
from the other, so that we can express fg and b/d as 
functions of these two parameters. In Figure 3, fg is 
plotted versus B — V and /i_B„. The gas fraction is 
larger for bluer colors B — V (panel a), which means 
that fg is larger for the MAHs whose present-day gas 
infall rate is still high. On the other hand, in panel 
b it is seen how a less concentrated disk (larger A) 
presents a higher fg than a disk with high S b„ (low 
A) . On the basis of this result is the infiuence of the 
disk gravitational compression upon the capability of 
gas to form stars (Firmani & Tutukov 1992, 1994). 
The disk surface density strongly influences the b/d 
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ratio. This is because the stellar surface density en- 
ters in the Toomre gravitational instability criterion 
which is used in our models to calculate the forma- 
tion of bulges. The larger the central SB (stellar 
density), the larger is the b/d ratio (Fig. 4b). The 
mass and MAH introduce a dispersion in this cor- 
relation because they influence the other quantities 
which appear in the Toomre criterion. In Figure 4a 
is seen how marginally does b/d depend upon B — V. 
The observational data show that b/d correlates bet- 
ter with fiBa than with B — V (e.g., de Jong 199()a). 

In Figures 3 and 4 are also plotted the observa- 
tional data taken from a compilation presented in 
McGaugh & de Block (1997) where low SB galax- 
ies are included. The B — V color indexes were 
not corrected for the internal (inclination) galaxy ex- 
tinction. We have applied this correction according 
to the formula given in the RC3 catalog (de Vau- 
couleurs et al. 1991). The b/d ratios for the low SB 
galaxies presented in McGaugh & de Block (1997) 
were not estimated, so that in the panels where this 
ratio is plotted the low SB galaxies are not consid- 
ered. In general, the model properties fall rather well 
within the observational ranges and are in agreement 
with the empirical correlations (compare also Table 1 
with the correlation matrix presented in McGaugh & 
de Block 1997). It is encouraging that the b/d ratios 
predicted by the models using the simple gravita- 
tional instability criterion (secular bulge formation) 
are in agreement with those inferred from observa- 
tions (de Jong 1996a, c). Note that de Jong (1996c) 
used an exponential profile in his two dimensional 
bulge-to-disk decomposition procedure instead of the 
de Vaucouleour's profile for the bulges, arguing that 
such a profile fits better the observations (see also 
Andrekakis & Sanders 1994; Andrekakis, Peletier, & 
Balcells 1995). That is why the b/d ratios obtained 
by him are smaller than those given by previous 
b/d decompositions (e.g., Simien & de Vaucouleurs 
1986). The bulges of galaxies with b/d ratios larger 
than those our models predict can be the product 
of the collapse and/or merger mechanisms of bulge 
formation which we do not treat. 

The most serious inconsistency between theory 
and observation in Figures 3 and 4 is probably that 
the color index of the models are bluer than those 
of the observed galaxies. Part of this inconsistency 
might be due to the imprecision of our population 
synthesis model (see §2). Nevertheless, other ex- 
planations are also possible. The statistical range 
(~ 80%) of the MAHs calculated here for the Gaus- 
sian fluctuations leads to disks with B — V between 
approximately 0.4 and 0.7. Models with colors red- 



der than 0.7 mag may be produced, although only for 
extreme cases. We have found that the color index 
becomes very sensitive to the MAH when this cor- 
responds to early active MAHs: for some extreme 
cases the B — V color may be as red as ~ 0.95 mag. 
Thus, some models can easily attain colors redder 
than 0.7 mag; of course the frequency of such mod- 
els will be low. In Figures 3 and 4, with dashed lines 
we show the range of values of the different proper- 
ties of the 5x10^^ Mq models for the three A's, when 
the statistical range in the MAHs is symmetrically 
extended to 94% (symbols consider only 80%). As 
is seen, the models can be as red as some observed 
galaxies are. 

The addition of some extra phenomena in our 
models may help to obtain galaxies redder than those 
presented in Figures 3 and 4. For example, the inter- 
nal extinction in combination with the metallicity- 
luminosity relation can introduce an important ef- 
fect of reddening, particularly for the most massive 
galaxies (see §5.2). Using the results presented in 
Wang & Heckman (1996) (hereafter WH96), and the 
Galactic extinction curve for Ry = 3.1 (Cardelli, 
Clayton, & Mathis 1989), we have reddened the 
models corresponding to the average MAH and A = 
0.05, represented in Figures 3 and 4 with black filled 
circles. The error bars account for the range of pa- 
rameter values given in WH96 (see Appendix). 

The influence of environment on the galaxy evo- 
lution might also help to produce models redder than 
0.7 mag. In dense environments, the gentle mass ag- 
gregation is probably truncated early. Experiments 
show that if gas accretion is truncated in the models 
at 6 Gyrs and 4 Gyrs, then B — V roughly increases 
by 0.08 and 0.15 mag, respectively. On the other 
hand, the early interactions in the dense environ- 
ments can induce non stationary SF which produces 
a fast gas consumption into stars. In a recent work, 
Marquez & Moles (1999) studied a sample of very 
isolated galaxies and they found that the corrected 
total B — V color index of the galaxies lies between 
0.35 and 0.85 magnitudes with an average value of 
0.48 magnitudes, i.e. the very isolated galaxies in- 
deed tend to be blue. 

The observational sample presented in Figures 3 
and 4 is statistically incomplete and plagued of ob- 
servational uncertainties. This and possible effects 
of extinction could account for the large scatters in 
the empirical correlations among the intensive prop- 
erties with respect to those of the models. Neverthe- 
less, it is probable that the intrinsic scatters are, in 
any case, larger than those found in our simulations. 
This is because the scenario proposed here takes into 
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account only the most relevant ingredients of galaxy 
formation and evolution, omiting certain phenomena 
which are not dominant, but are observed in some 
galaxies. For example, the SF prescription used in 
the models is basically a stationary process, while in 
real galaxies SF sometimes may appears in bursting 
modes. This fact introduces a stochastic component 
in the photometric features, particularly for low mass 
galaxies (Firmani & Tutukov 1994). 

5. CORRELATIONS BETWEEN GLOBAL 
PROPERTIES 

5.1. The Tully- Fisher relation in the infrared bands 

and its scatter 

The galaxy luminosity in the infrared bands is a 
reliable tracer of the disk stellar mass (e.g.. Pierce & 
Tully 1992; Gavazzi 1993; Gavazzi, Pierini, & BoscUi 
1996). Thus, the TFR in the infrared bands may be 
translated to a relation between the disk stellar mass 
Ms and the maximum rotation velocity Vm of disk 
galaxies. In AFH98 it was already shown that this 
relation for the ag ~ 0.6 SCDM model is in good 
agreement with that derived from most of observa- 
tions. AFH98 assumed that the total disk mass is 
directly proportional to Mg. Since here we calculate 
the SF, this assumption is relaxed. The Alg — Vm 
relation we obtain is actually very similar to that 
estimated in AFH98. However, the scatter in the re- 
lation now is smaller than in AFH98. This is because 
Mg depends upon disk surface density. The efficiency 
of SF is larger in disks with larger surface densi- 
ties. Thus, for a fixed disk total mass (stars+gas), 
the model galaxies with higher surface densities have 
larger Vm than models with low surface densities, 
but they also have larger stellar masses Af,. In the 
Ms — Vm diagram, the models in the high velocity 
side shift to high values of Ms, while the models in 
the low velocity side shift to low values of Ms. This 
compensating effect is responsible for a scatter in the 
Ms — Vm relation smaller than that we would expect 
if one assumes that Ms is directly proportional to the 
total disk mass. This effect also produces that high 
and low surface density (brightness) models have ap- 
proximately the same Mg — Vn relation, i.e. the TFR 
of low and high models is nearly the same (see also 
FA). The observational data point out to the same 
fact (c.f., Zwaan et al. 1995; Verheijen 1997). 

In our models the rms scatter in the Ms — Vm 
relation (or TFR) is originated by the scatter in the 
structure of the DM halos (due to the scatter in the 
MAHs) , and by the dispersion of the spin parameter 
A. The contribution of the latter to the total scat- 
ter is diminished due to the effect just mentioned 



above. To estimate the scatter in the Ms — Vm rela- 
tion we should generate a catalog of models where, 
at least, the MAH and A are drawn through Monte 
Carlo simulations from their corresponding statisti- 
cal distributions. This was done in FA for a low 
density CDM model with cosmological constant. At 
galaxy scales the power spectrum of fluctuations of 
this model is similar to the SCDM model used here. 
In FA it was obtained that the quadratic contribu- 
tions to the scatter in the Ms — Vm relation due to 
the scatters in the MAHs and in A are approximately 
75% and 25%, respectively. Thus, the scatter in the 
MAH is the main source of the scatter in the Mg — Vm 
relation (or TFR). 

It is important to note that this scatter is corre- 
lated with some intensive properties. For example, 
in Figure 5 is seen that, for a given mass, Vm in- 
creases with the B — V color. This is also evident 
from Table 1 where we show the correlation coeffi- 
cients for Atf = Li/V™' which reflects the devia- 
tions from the TFR. This behaviour is easily under- 
stood from the point of view of the extended collapse 
scenario: galaxies formed through gentle MAHs will 
be less concentrated (smaller Vm) and with SF histo- 
ries more extended in time (bluer colors) than galax- 
ies formed through early active MAHs. The obser- 
vational data confirm this prediction of the models. 
In Figure 5, galaxies from a cross of the RC3 (de 
Vaucouleurs et al. 1991) and the Tully (Tully 1988) 
catalogs are also depicted. To estimate the behav- 
ior of Vm with the B — V color for a given mass 
(luminosity), the data were divided into three broad 
bins according to the S— luminosities presented in 
the mentioned catalogs. The dashed lines are lineal 
regression to each one of these bins. Although the 
estimate is very qualitative, we see that the observed 
trend is in agreement with the models' predictions. 

5.2. The color Tully-Fisher and color-magnitude 

relations 

The measured slope of the TFR is smaller as the 
color passband in which the luminosity is measured 
is bluer. This dependence of the TFR slope with 
the color band gives rise to the so called color TFR 
(Gavazzi 1993; Bothun et. al 1985): 

l^iocVl--'"^^ (4) 

where i and j indicate a given color band, and mj 
and TOj are the corresponding slopes of the TFR in 
these bands. For example, for the H and B bands, 
according to the empirical estimates reported in the 
literature, (mjj — ms) w 0.4 — 1.2 (e.g., Gavazzi 
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1993; Strauss & Willick 1995, and more references 
therein). As it was mentioned in the previous sub- 
section, AIs oc Lh- Thus, according to eq. (4) we 
have that 

Ms/Lb cx Vi, (5) 

where /3 = tuh — ruB ~ 0.4 — 1.2, i.e. the Ms/Lb 
ratio increases with Vm- This is an important re- 
sult which should be seriously considered in any the- 
ory of disk galaxy formation. We find at least three 
possible causes for this dependence of the mass-to- 
luminosity ratio on Vm' (i) the SF efficiency could 
depend on the gravitational potential of the system 
in such a way that in galaxies with larger Vm the 
gas is consumed into stars more efficiently than in 
galaxies with low Vm', (ii) if disk galaxies build up 
inside-out by a continuous process of gas accretion, 
then the dependence given by eq. (5) is expected 
in the case the rate of gas accretion at late epochs 
decreases more rapidly for larger galaxies than for 
smaller ones^ (iii) the color TFR or the dependence 
given by eq. (5) may be explained by the effect of 
internal dust absorption — which affects the luminos- 
ity in the blue bands but not in the infrared ones — , 
if the absorption increases with the galaxy mass (ro- 
tation velocity). 

The results of our models show that the slopes 
iriH and tub of the TFRs in the B— and 77— bands 
are approximately equal, or, what is the same, the 
Ms/ Lb ratio is constant with Vn- This is (i) be- 
cause, according to our SF mechanism, the SF ef- 
ficiency scarcely depends on mass (or Vm), and (ii) 
because we assume that the rate of gas accretion 
is determined only by the cosmological rate, avoid- 
ing possible intermediate processes like that men- 
tioned in the footnote. In fact, in our models, at 
late epochs the cosmological accretion rates — given 
by the MAHs — for larger systems decrease more 
slowly than for smaller systems (small halos collapse 
on average earlier than big halos, Navarro, et al. 
1996, AFH98). However, this dependence on mass 
(oposite to observations) is compensated by the fact 
that larger galaxies have on average slightly higher 
surface densities (see §3; Dalcanton et al. 1997; FA), 
and in our models, the SF efficiency is larger for 
higher surface densities. It is interesting to note 
that the same argument was used to explain why 
the shape of the rotation curves is nearly indepen- 
dent of mass in spite of that less massive CDM halos 
are more concentrated than the more massive ones 

^Such a situation is possible if, for example, in large sys- 
tems a fraction of the accreting gas is trapped and efficiently 
transformed into stars within the potential wells of relative 
large substructures (satellites) around the central disk galaxy. 



(§3). 

We propose that the dependence given by eq. 
(5) is basically produced by dust absorption. It is 
well known that the metallicity abundance of spiral 
galaxies correlates with their luminosities, Z oc L%, 
where e w 0.3 — 0.5 (see references in Roberts & 
Haynes 1994). Since dust particles form from heavy 
elements, this relation suggests that more luminous 
galaxies have higher dust abundance than the less 
luminous galaxies. Direct galaxy images also show 
that the dust abundances increase with the galaxy 
luminosity (van den Bergh & Pierce 1990). Indeed, 
WII96 based on studies of the UV and FIR fluxes of 
a sample of normal late-type galaxies, have found a 
clear correlation between the dust opacity and the 
_B— band luminosity of galaxies. 

In the Appendix, from the results of WII96 and 
using the uniform slab model, we calculate the ex- 
tinction in the i3— band, Ab, as a function of Lb 
(eq. A2). Applying this luminosity-dependent ex- 
tinction to our models we obtain that the slope of 
the TFR in the B— band becomes flatter than in the 
original models. Figure 6 shows that this relation in 
the range logLs / Lbq ~ 9.5 — 11.5 is approximated 
by a line with slope ~2.7, as many observations sug- 
gest. Note that the dependence of extinction on lu- 
minosity does not only produce a change of slope, 
but also some nonlinearity, particularly at the bright 
end of Figure 6. A similar result was previously re- 
ported by Giovanelli et al. (1995). If the samples 
used by different observers cover different luminos- 
ity ranges, this could explain why the reported slopes 
of the B-band TFRs are so disparated. Recently, 
Kudrya et al. (1997) have presented the B— band 
TFR for a large sample of galaxies in a large range 
of luminosities; from their Figure 6 it is clearly seen 
how the slope of this relation tends to be steeper 
for the less luminous galaxies (see also e.g.. Pierce 
& TuUy 1992). In Figure 6 are also depicted the 
predicted TFRs after applying extinctions calculated 
with the maximal and minimal optical depths given 
in WII96 (see Appendix). The more realistic "sand- 
wich" model (Disney, Davies, & Philipps 1989) was 
also considered (dotted line, see Appendix). 

Since the early eighties a correlation between the 
galaxy total color and its magnitude was reported for 
spiral galaxies (Vishnavatan 1981; Wyse 1982; TuUy, 
Mould, & Aaranson 1982). The more luminous is 
the galaxy the redder tends to be its color. As a 
matter of fact, this relation is intimately related to 
the color TFR mentioned above. For instance, by 
dividing the B— band Lb = CbV///^ and i/— band 
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Lh = CrV^" TFR one to another, one obtains: 
= 2.5(^-l)logLB+2.51og f 

(6) 

For the same range of values for ms and niH used 
in eq. (5), we have that {B — H) oc a log Lb with 
a ?a 0.4 — 1.3. This is in rough agreement with 
the slopes reported for the color-magnitude relation. 
Thus, according to the above discussion, dust ab- 
sorption may be responsible for this relation. Assum- 
ing that in the if— band dust absoption is negligible, 
while in the B— band absorption depends upon lumi- 
nosity as WH96 found (cq. Al), the galaxy {B — H) 
color will redden proportionally to the extinction Ab 
estimated in the Appendix. As a first approxima- 
tion, using the linear fitting to (eq. A3), one 
obtains that 

(B-ff) oc 0.42 log Lb- (7) 

WH96 suggested that the dependence of inter- 
nal absorption upon luminosity (cq. Al) may be 
explained by the observed increase of metallicity 
and/or surface density with the luminosity. Less lu- 
minous (massive) galaxies can have lower metallic- 
ity abundances than more luminous galaxies because 
the dust expelling from the galaxy is more effective in 
the former than in the latter. Shustov, Wiebe, & Tu- 
tukov (1997), using models of disk galaxy evolution 
that take into account the process of dust expelling 
by radiation pressure, show that dust loss is much 
more effective in low-mass galaxies than in massive 
ones. They obtained a dependence of iron abundance 
on galaxy luminosity near to the observed, and found 
that the other properties of the galaxy models do not 
significantly change with dust loss. Hydro dynamical 
simulations of starburst-driven mass ejection indeed 
have shown that the efficiency of metal ejection is 
much larger than that of the gas and this efficiency of 
metal ejection decreases with mass (Mac Low & Fer- 
rara 1999). Regarding the disk surface density, our 
models indeed show that more massive galax;ies have 
typically slightly higher surface density disks. Al- 
though this dependence is not too pronounced as the 
one Dalcanton et al. (1997) obtained, it could also 
contributes to give rise to the absorption-luminosity 
relation (eq. Al). 

Finally, we note that metallicity, in addition to 
the absorption effect, also influences the spectropho- 
tometric evolution. For example the spectrophoto- 
metric models of Bressan, Chiosi, & Fagotto (1994) 
applied to single stellar populations with different 



initial metallicities show that at ~ 12 Gyr the dif- 
ferences in the V — K and B — V colors are ~0.88 
mag and ~0.23 mag, respectively for a factor of 20 of 
variance in metallicity. These models roughly agree 
with the observed color-magnitude relation of ellipti- 
cal galaxies, and in accordance with Kodama & Ari- 
moto (1997), this relation is consequence of a metal- 
licity effect, instead of an age effect. 

5.3. The luminosity-radius relation 

The mass-radius relation predicted by the models 
(the average values in the MAH and A were used) is: 

Ms oc (8) 

where Mg is the disk stellar mass, and Rr is the 
Holmberg radius. The radius scales with Lb as R?jf- 
The scatter in these relations is large and it corre- 
lates with the SB (AF99). The Holmberg and scale 
radii do not correlate with any intensive galaxy prop- 
erty (see Table 1), suggesting that the evolution of 
disk galaxies and the Hubble sequence are size in- 
dependent (de Jong 1996a; McGaugh & de Block 
1997). 

6. THE FUNDAMENTAL PHYSICAL FACTORS 
OF DISK GALAXIES AND THE HUBBLE 
SEQUENCE 

One of the aims of this paper is to investigate 
which physical factors determine the properties of 
present-day disk galaxies. According to the extended 
collapse scenario, these fundamental factors are: the 
virial present-day mass My of the system, the MAH, 
and the angular m,omentum expressed through the 
spin parameter A. These factors and their statisti- 
cal distributions are related to the initial cosmologi- 
cal conditions and they are able to produce most of 
the disk galaxy properties and their correlations. As 
a first approximation, here we assumeed indepen- 
dence between A and the MAH. Analytical studies 
(c.f., Hoffman 1986, 1988; Heavens & Peacock 1988) 
suggest that A is almost independent from the fluc- 
tuation peak height (related to the MAH) for CDM 
power spectra. We also assumed independence be- 
tween A and My (see §2.3). Concerning the depen- 
dence of MAH upon My, the less massive galaxies 
undergo a faster early collapse than the more mas- 
sive galaxies for the CDM power spectra (AFH98). 

The correlation coefficients from a principal anal- 
ysis of these three factors with the model galaxy 
properties are given in Table 2. The MAH was quan- 
tified through the 7 parameter, where 

^ log{Mnom/{Mno„,/2)) 
^ /05(i(M„om)A(M„oxn/2))' ^ ' 
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The MAH, which partiaUy drives the SF history, 
strongly influences the B — V color index, and mod- 
erately influences the gas fraction fg. The range of 
B — V colors that the models span is mainly asso- 
ciated with the statistical dispersion in the MAHs. 
This dispersion is also reflected in the TFR scatter 
(the coefficient A^p correlates with 7) , and as it was 
shown above (see Figure 5), observations confirm a 
correlation between B—V and the maximum circular 
velocity for a given luminosity (mass). 

The A parameter strongly infiuences the SB, the 
b/d ratio, and fg. For a given My and MAH, A deter- 
mines the degree of concentration (SB) of the disks. 
According to the SF mechanism used in our models, 
the surface density influences the SF efficiency and 
therefore, the fg. Furthermore, as it was mentioned 
in § 4, the stellar Toomre parameter Qs is small for 
small stellar surface densities. Therefore, the b/d 
ratio is larger for systems with small A's. 

Mass strongly correlates with luminosity and the 
Holmberg radius and slightly influences the b/d ra- 
tio and the SB (both in the same direction as sug- 
gested by observations). As it was pointed out in 
§3.2, B — V and /is^ are the two parameters on which 
the other intensive properties depend, i.e. the inten- 
sive properties of disk galaxies may be described in a 
biparametrical sequence, whose origin deals with two 
of the fundamental physical factors of galaxies, the 
MAH and A, respectively. Moreover, we have found 
that these properties correlate between them follow- 
ing the same trends that observed galaxies present 
across the Hubble sequence. Most of the intensive 
properties are almost independent from the third 
fundamental factor, the mass (luminosity). We have 
found that the mass-dependent dust absorption may 
be responsible for the color TF and color-magnitude 
relations (§5.2). 

The morphological Hubble classification has been 
a useful guide to study the observational proper- 
ties and correlations of galaxies. Nevertheless, most 
of the classification criteria deal with morphological 
characteristics that probably are transient phenom- 
ena related to more fundamental galaxy properties. 
The main classification criteria for spirals are the 
pitch angle and the strength of the spiral arms, and 
the b/d ratio. Despite the big effort done, the prob- 
lem of the origin and maintenance of arms in disk 
galaxies is still not well understood. Bertin & Romeo 
(1987), Bertin et al. (1989a,b) and other authors 
(see for references Combes 1993) coincide in point- 
ing out that gas, dynamically speaking, is crucial for 
the excitation and maintenance of spiral structures. 
According to Bertin & Lin (1996), the gas fraction 



mainly determines the sequence of types a, b, and 
c, i.e. the pitch angle of arms. As is seen in Table 
1, the gas fractions in our models correlate with the 
other secondary indicators of the Hubble sequence 
which go across the a, b, and c sequence, such as 
b/d, and B — V, suggesting that the origin of 
the Hubble morphological types is closely related to 
the galaxy formation and evolution processes of the 
extended collapse scenario. 

Several observers have pointed out that, for a 
given luminosity, galaxies with larger maximum ro- 
tation velocities (the TFR scatter) are of earlier 
types than those with smaller velocities (Roberts 
1978; Rubin, Thonnard, & Ford 1980, Rubin et 
al. 1985; Giraud 1986, 1987; Krann-Korteweg, 
Cameron, & Tamman 1988; Giovanelli et al. 1997). 
The TFR, Li = Atf^C"' where i is a given band, 
predicted by the models has a dispersion that, at a 
first approximation, we express only through vari- 
ations in the coefficient Atf- As has been com- 
mented in § 3.2 (see also Fig. 5), the models present 
correlations between j4tf and some secondary Hub- 
ble type indicators as B — V, fg, and b/d (Table 
1). Thus, within the framework of the extended col- 
lapse scenario, the empirical correlation found be- 
tween galaxy type and rotation velocity for a given 
luminosity finds a natural explanation. 

7. DISCUSSION 

7.1. Problems of the hierarchical formation scenario 

We have studied disk galaxy formation and evo- 
lution in the cosmological context taking into ac- 
count several key factors and processes of this phe- 
nomenon. Our models represent an extreme case of 
galaxy formation where we have assumed: (i) gentle 
and spherical symmetric gas infall at the center of hi- 
erarchically growing CDM halos (no mergers), (ii) no 
gas reheated and expelled back into the halo or lost 
from the system, and (iii) stationary self-regulated 
SF within the disk. We have seen that at the level 
of this simple picture many properties and correla- 
tions of the observed disk galaxies were predicted. 
However, the models are in potential conflict with 
some observational pieces of evidence. As we shall 
see, these conflicts may be partially solved invok- 
ing some extra ingredients in our scenario, but it is 
also possible that they are pointing out to serious 
inconsistences of the hierarchical clustering scenario 
in general. 

In an accompanying paper (Firmani & Avila- 
Reese 2000b) we study the evolution of the models 
presented here. We find that the size and SB evolu- 
tion of the disks is more dramatic than that one the 
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deep field observations seem to suggest. This prob- 
lem, toghether with the too steep negative radial 
color gradients our models predict (§3), constitute 
serious difRcultics for the extended collapse scenario. 
The color gradients we obtain might be even steeper 
if the galaxy's internal extinction significantly de- 
pends on surface density. The model disks also tend 
to be somewhat less dense and more extended than 
real galaxy disks, and the gas fractions of the mod- 
els are slightly larger than those measured in galaxies 
(see Figs. 3 and 4). On the other hand, the integral 
B — V colors of observed red galaxies seem to be red- 
der than the red models (Figs. 3 and 4). In §4 we 
have extensively discussed the solutions to the last 
conflict. 

If some transfer of angular momentum from 
baryon matter to DM is allowed, then the disks re- 
sult smaller and more concentrated. If this transfer 
of angular momentum tends to be more important 
at latter epochs and/or if a fraction of the late, high 
angular momentum accreting gas is trapped within 
substructures and does not fall to the disk, then the 
color gradient will be smaller, the total color index 
will be redder, and the size evolution will be less pro- 
nounced than we have obtained. It is also possible 
that during the gas collapse, the original gas angular 
momentum is redistributed due to several dynamical 
and hydrodynamical processes related to the devia- 
tions from the homogeneity and sphericity. Exper- 
iments where a moderate "flattening" of the radial 
angular momentum distribution of the infalling gas 
was introduced, show us that the disk radial color 
index gradients significantly decrease and the disk 
sizes change more slowly in time with respect to the 
original models. 

It is interesting to note that the above mentioned 
conflicts of the extended collapse scenario tend to 
be opossed to those of the merger scenario. The 
merging seem to cause serious difficulties for con- 
serving the gas initial angular momentum, with the 
consequence that disks result too small. This prob- 
lem, called the "angular momentum catastrophe" , 
has been observed in the N-body-|-hydrodynamical 
numerical simulations without SF (e.g., Navarro & 
Steinmetz 1997). When SF is included, the merging 
of the large stellar subunits as the system evolves, 
also produces loss of angular momentum with the 
consequent formation of a elliptical-like object in- 
stead of a disk (c.f., Katz 1992). It is probable 
that reality lies in between the extended collapse and 
merger scenarios. More observational and theoreti- 
cal studies on the formation and evolution of disks, 
their SF histories, and the radial distribution of color 



indexes and extinction are needed in order to arrive 
at definitive conclusions. 

While some physical ingredients not considered 
in our scenario actually might work in the direction 
of improving the predictions, we do not discard the 
possibility that in general the hierarchical formation 
picture is in conflict with reality. On the other hand, 
as already has been pointed out, if the merging ac- 
tivity in the hierarchical formation picture based on 
the Gaussian inflationary CDM models is actually 
too high, then it is difficult to figure out the exis- 
tence of a large number of galaxy disks in the frame- 
work of this picture (e.g., Toth & Ostriker 1992). To 
the problems related with the formation of disks in 
the standard hierarchical formation picture, we add 
other potential conflict: the inner structure of the 
CDM halos. This is the subject of the next subsec- 
tion. 

7.2. Effects of a shallow core in the DM halo 

The rotation curves of dwarf and low SB galax- 
ies suggest that the inner density profiles of DM ha- 
los are shallower than (Moore 1994; Flores & 
Primack 1994; Burkert 1995; de Blok & McGaugh 
1997). Our results and the results of high-resolution 
cosmological N-body simulations (e.g., Navarro et al. 
1996,1997; Fukushige & Makino 1997; Moore et al. 
1998, 1999; Jing 1999; Jing & Suto 1999) show that 
for the CDM models the halo inner density profiles 
are not shallower than r~^. The rotation curve de- 
composition of the galaxies formed within the CDM 
halos also seem to be in conflict with the decomposi- 
tions inferred from observations for high SB galaxies 
(AFH98; FA; but see Navarro 1998). 

In order to explore the effects of a shallow core on 
the properties of the model disk galaxies, we intro- 
duce an artificial core in our evolving CDM halos. 
Our approach is purely phenomenological: the in- 
ner density profile of the evolving halos is fiattened 
in such a way that the final rotation curves of low 
SB model galaxies (those with large A) resemble the 
observed rotation curves of low SB (see for details 
FA). The inner halo density profile obtained for the 
SCDM cTg = 0.6 model treated here is very similar 
to that obtained for the ACDM model in FA. 

As in FA, we find that the main effect of the 
shallow core is to raise the contribution of the disk 
component with respect to the halo component in 
the rotation curve decomposition. This makes the 
rotation curve decompositions of the model galaxies 
more realistic. The intensive properties are almost 
the same for the models with and without a shallow 
core. The Vm of the models with a shallow core 



14 



AVILA-REESE & FIRMANI 



arc slightly smaller than for the original ones. This 
shifts our predicted /— and ff— band TFRs to lower 
velocities thus giving a better agreement with the 
observational data (sec also FA) . For a fixed rotation 
curve shape, the central surface densitites (SB) are 
higher for the models with core. This is because in 
order to have the same rotation curve shape, a model 
with a shallow core should have a more concentrated 
disk (smaller A) that would compensate the smaller 
halo contribution to the total rotation curve with 
respect to the model with a more dense halo core. 

8. SUMMARY AND CONCLUSIONS 

We studied the formation and evolution of in- 
dividual disk galaxies within the framework of an 
inflationary Gaussian SCDM model normalized to 
CTg = 0.6. We built-up disks in centrifugal equilib- 
rium within hierarchically growing DM halos; the 
disks form inside-out with gas (no mergers) accreted 
with a rate proportional to the cosmological infall 
rate given by the mass aggregation history (MAH), 
and under the assumption A = const in time and de- 
tailed angular momentum conservation during the 
gas collapse. The SF in the disk is produced by 
global gravitational instabilities and it self-regulates 
by an energy balance in the turbulent ISM. We cal- 
culated the secular formation of a bulge using a grav- 
itational instability criterion for the stellar disk. Our 
main conclusions are: 

1) . The main properties of present-day high and 
low SB disk galaxies and their correlations may prin- 
cipally be explained by the combination of three fun- 
damental physical factors and their statistical distri- 
butions. These factors are the halo virial mass M^, 
the MAH, and the spin parameter A, and they are 
related to the cosmological initial conditions. 

2) . The SF efficiency of the models does not sus- 
tantially depend upon the mass or Vm of the galaxy. 
Therefore, the integral color index B—V or the mass- 
to-luminosity ratio M/ Lb do not depend upon mass. 
We showed that the empirical correlation of dust ex- 
tinction upon luminosity reported by WH96 explains 
the observed color-magnitud and color TFRs with- 
out the necessity of introducing a mass-dependent 
SF efficiency. 

3) . Several global and local properties of present- 
day disk galaxies and their correlations are a natural 
consequence of the extended collapse scenario. We 
remark the following results: 

- The redder and more concentrated is the model 
disk, the smaller is its gas fraction and the larger is 
its b/d ratio. The values of these intensive proper- 
ties are in rough agreement with the observational 



data. The intensive properties of the models are de- 
scribed by a biparametrical sequence, where the pa- 
rameters are the color index B — V and the central 
SB, fJ-Bo- These parameters arc mainly determined 
by the MAH and by A, respectively. Thus, according 
to our models, the Hubble sequence is biparametrical 
as some observations suggest. 

-The model results agree with the near-infrared 
TFRs, suggesting that the origin of the TFR is cos- 
mological. This relation is almost independent from 
the SB and its scatter is mainly dominated by the 
dispersion in the MAHs (or the dispersion in the DM 
halo structures). For a given luminosity, the larger 
the Vm is (TFR deviation), the redder is the galaxy 
and the smaller is its gas fraction fg. This is because, 
on one hand the MAH influences on the galaxy color 
index and fg, and on the other hand, it is the dis- 
persion in the MAHs that mainly produces the de- 
viations from the TFR. 

-The stellar surface density and S— brightness 
proflles are nearly exponential. This is a result of the 
inside-out disk formation assuming detailed angu- 
lar momentum conservation and A constant in time. 
The SB of the models mainly depends on A. 

-Also result of the inside-out disk formation, 
the radial B — V gradients along the disks are neg- 
ative: the outer disk regions are younger than the 
inner ones. 

-The shape of the total rotation curves depends 
on the the disk mass fraction fd, A and the MAH. 
For reasonable values of fd, the approximate flatness 
of the rotation curves presented by most of the mod- 
els is due to a "conspiration" among the CDM halo 
and the baryon disk. The rotation curve decomposi- 
tions, however, show excesive dominance of the halo 
component over the disk component with respect to 
observational estimates. 

The main shortcomings of the models presented 
here are (i) the probably excesive radial color index 
gradients and disk size evolution, and (ii) the DM 
dominion in the rotation curve decompositon. The 
inclusion in the models of other, less important ingre- 
dients and processes as merging, angular momentum 
transfer, and non-stationary SF work in the direction 
of overcoming the former shortcomings. The latter 
problem is solved if the inner density profile of the 
CDM is shallower than predicted. It is probable that 
these problems are in general pointing to a more se- 
rious difficulty of the hierarchical formation scenario 
based on the Gaussian CDM models. 
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A. CALCULATION OF THE 
LUMINOSITY-DEPENDENT B-BAND 
EXCTINCTION 

Wang & Heckman (1996), based on studies of the 
far UV and FIR fluxes of a sample of normal late 
type galaxies, have concluded that the dust opac- 
ity increases with the luminosity of the young stellar 
population. They find that a power law relation be- 
tween the optical depth and the UV luminosity ex- 
plains this observational dependence. This relation, 
referred to the B band, is: 



TB = TB,. 



Lf 



■IB,* 



(10) 



where the best fits to observations are for Lb.* = 
1.3 X 101"Lbq , TB,* = 0.8 ± 0.3, and (3 = 0.5 ± 0.2. 
According to the uniform slab model WH96 have 
used, the extinction in magnitudes may be expressed 
as Ab = —2.5 log ^ ^~g^P^~'^^) j ^ and in the range 
IQ-^^Lfl^, using (Al) with the central values, 



10^ 



is well approximated by 



Ab « 0.38+0.42 log ( -J^^] +0.14 (log ( 



IQIOLb, 



lOiOLi 



In the range of 10^ — 10^^Lb„ the linear fitting 



(11) 



Ab « 0.43 + 0.42 log 



(12) 



is also a good approximation. 

To calculate the extinction for the two extreme 
cases of maximal and minimal optical depths, we use 
tb,* = 1.1 and (3 = 0.7, for the former, and tb,* = 
0.5 and (3 = 0.3 for the latter. To calculate the 
extinction with the more realistic "sandwich" model 
(Disney, Davies, & Phihpps 1989), we use the fiducial 
value of (i given by WH96, and tb,* and the ratio of 
the height scale of dust to young stars, C„ are fixed 
to the values suggested by Bosselli & Gavazzi (1994) 
for the optically thin case in the H band (tb,* =1.33 
and C =0.74). 
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Figure 2. Rotation curves for the same models of 
Figure 1 (5xlO"Af0), and for a 5x10^" Mq (bot- 
tom curve) and SxlO^^M© (top curve) galaxy corre- 
sponding only to the average MAH, A = 0.050. Radii 
were scaled to the optical (Holmberg) radii of each 
model. 



Figure 1. The B-surfacc brightness (a) and the B — 
V color index (b) profiles of a SxlO^^M© galaxy. 
The average MAH cases for the spin parameters A = 
0.035 (dashed line), A = 0.050 (solid line), and A = 
0.100 (point-dashed line) are represented with the 
thick lines. For the early active MAH only the model 
with A = 0.035 (thin dashed line) is plotted, while 
for the extended MAH, only the A = 0.100 case (thin 
point-dashed line) is plotted. All the models were 
calculated for a erg = 0.6 SCDM model. 
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Figure 3. The gas fraction fg vs. the integral B — V color index (a), and vs. the central B-SB (b) for 
models and observations. The gray, black, and white filled symbols, correspond to models with A = 0.035, 
A = 0.050, and A = 0.100 respectively. Squares are for the early active MAH, circles for the average MAH, 
and triangles for the extended MAH. Three masses (dark+baryon), 5xlO^°Mo, 5xlO^^M0, and 5xlO^^M0 
are considered (the larger the mass, the smaller is the gas fraction). The dashed lines connect the models of 
constant mass for SxlO'^-'^M©, and extend the statistical range of MAHs to 94% (symbols consider only 80% 
of the MAHs). The three small black filled circles are the same models corresponding to the big black filled 
circles but reddened according to the dust absorption-luminosity dependence given in WH96 (see text). The 
error bars correspond to the range of values given in WH96. Small empty circles arc the observational data 
collected by McGaugh & de Blok (1997) and corrected for inclination extinction. LSB galaxies are included. 
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Figure 4. The b/d ratio vs the integral 5 — F color index (a), and the central B-SB jiBo (b) for models and 
observations. The same symbol and line codes of Figure 3 are used. The bulge-to-disk ratios were taken from 
the iC— band two-dimensional decompositions carried out by de Jong (1996b). LSB galaxies and a few normal 
galaxies shown in Figure 4 are absent in this Figure. 
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Figure 5. The maximum rotation velocity, vs. B — V for models and observations. The same symbol codes 
of Figure 3 are used. The observational data (small symbols) were taken from a cross of the RC3 and the 
TuUy (1988) catalogs (sec text). The small triangles, diamonds, and inverse triangles correspond to galaxies 
with luminosities in B band within the 10® - SxIO^Lbq, 3x10^ - 3 x lO^^isg, and 3x10^° - 2 x 10"Lbq 
ranges, respectively. The dashed lines are linear regressions to the observational data corresponding to these 
ranges. Note how the maximum velocity of models and observations for a given mass (or range of luminosities) 
correlates with the B — V color. 
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Figure 6. The predicted B-band TFR for the as = 0.6 SCDM model (thick sohd hne). The slope of this 
relation is 3.5. The other black lines show how the intrinsical TFR is transformed if the -B-luminosities are 
dimishcd by dust absorption according to the observational dependence of optical depth of dust on luminosity 
given in WH96. While the dashed line corresponds to the fiducial optical depth, the point-dashed and two 
point-dashed lines are for two extreme cases of maximal and minimal optical depths (see text). The point line 
was obtained using a sandwich model with a ratio of height scale of dust to young stars of 0.74 (see appendix 
for references). In the range logLs/ ~ 9.5 — 11.5 the dashed line is well approximated by a line with slope 
~ 2.7. The dotted and dashed gray curves are the empirical TFRs given in Kudrya et al. (1997) and Pierce & 
TuUy (1992), respectively. We have assumed W^^"^"^ = 2 x Knax- These TFRs have a steeper slope for velocities 
lower than Vmax ~ 100 km/s (see Figure 6 of Kudrya et al. and Fig. 3 of Pierce & TuUy). 
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Table 1 

CORRELATION MATRIX OF THE GLOBAL 
PROPERTIES 
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-0.07 
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-0.05 
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-0.82 
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-0.87 


0.27 


-0.45 
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0.68 
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-0.38 


0.54 


0.10 


-0.51 


Mb 










-0.94 


-0.95 


-0.20 














0.84 


-0.04 


log hd 














0.34 



(a) j4tf = Li/V^' reflects the deviations of the mod- 
els from the TFR 



Table 2 

CORRELATION MATRIX OF THE GLOBAL 
PROPERTIES AND THE FUNDAMENTAL 
PARAMETERS 
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Mb 




log hd 


Atf 


logM^"^ 


-0.44 


-0.24 


-0.17 


0.27 


-0.9 


0.92 


0.98 


0.25 




0.01 


-0.95 


0.48 


-0.12 


-0.24 


0.05 


0.29 


0.82 


A(c) 


0.86 


-0.07 


0.81 


-0.88 


0.13 


-0.23 


0.51 


0.34 



(a) nominal present-day mass 

(b) MAH parameter (see eq. 9) 

(c) spin parameter 



